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' Abstract. Wc prove a compact embedding theorem in a class of spaces of piecewise 

' functions subordinated to a class of shape regular, but not necessarily quasi-uniform trian- 

. gulations of a polygonal domain. This result generalizes the Rellich-Kondrachov theorem. 

["t I , It is used to prove generalizations to piecewise functions of nonstandard Poincare-Friedrichs 

' inequalities. It can be used to prove Korn inequalities for piecewise functions associated with 
elastic shells. 
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1. Introduction 

^ I Let f2 be a bounded polygonal domain in M^. Let Th be a shape regular, but not nec- 

I essarily quasi-uniform, triangulation on Q. We also use Th to denote the set of all (open) 
I triangular elements of the partition. Let Hi be the space of piecewise functions subordi- 
CS| ! nated to the triangulation Th- A function in Hj^ is independently defined on every element 
^ ! T E Th on which it belongs to H^{t). This kind of space of piecewise functions arises in 
^ . analysis of discontinuous finite element methods. It is desirable to generalize the Sobolev 
■ space theory to such spaces. Poincare-Friedrichs type inequalities have been proposed and 
^ ■ proved in several forms in the literature [21 HJ [8]. In this paper, we prove a compact embed- 
ding theorem that generalizes the Rellich-Kondrachov theorem. Such compact embedding 
theorem can be used, together with a modified compactness argument, to prove the afore- 
mentioned Poincare-Friedrichs type inequality for piecewise functions in a general form, 
including those in nonstandard forms [7]. It seems necessary to generalizing Korn's inequal- 
ities on curved surfaces [3] to spaces of piecewise functions, which is useful in analyzing 
discontinuous Galerkin finite element methods for models of elastic shells. 

Shape regularity of triangulations is a crucial notion in this paper. It is worthwhile to 
recall its definition here. Considering a triangle, we let r and R be the radii of its largest 
inscribed circle and smallest circumscribed circle, respectively. Then the ratio R/r is called 
its shape regularity constant, or simply shape regularity. For a triangulation, the maximum 
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of shape regularity constants of all its triangles is called its shape regularity constant [6], 
denoted by /C. We will need to consider a (infinite) family of triangulations. For a family, 
the shape regularity constant /C is the supreme of all the shape regularity constants of its 
triangulations. For a given triangulation Th, we use Qh to denote the union of all the open 
triangular elements, and use £^ denote the set of all interior (open) edges and all boundary 
edges. A function u in Hj^ is certainly in L^(f2) (in which the function value on has no 
significance). On an edge e G S^, a function u has two different traces from the two elements 
sharing e. We use |u] to denoted the difference of the two traces, which is the jump of u 
over e. We furnish the space Hj^ with the norm 



\u\ 



H, 



iVnl 



1/2 



Here |e| is the length of the edge e. The integrals are taken with respect to Lebesgue measures 
of the integration domains. Let {Th^} be a family of triangulations with a certain regularity 
constant /C. Let Ui G be uniformly bounded sequence of piecewise functions, i.e., 
there is a constant C such that < C for all i. Then there is a subsequence {ujj.} 

that is convergent in L^(f2). It is this compact embedding theorem on families of spaces that 
is needed to prove, for example, a Poincare-Friedrichs inequality that there is a constant 
C that is dependent on VL and dependent on the triangulation Th only through its shape 
regularity but otherwise independent of Th such that 

1/2 



|«||L2(n) < C 



\Vu\ 



^ e 



u 



Vug Hl{n). 



Such constant independence of triangulations is fundamental in numerical analysis. 

The remaining of the paper is organized as follows. In Section |2], we present a trace 
theorem for piecewise functions. This result will be used several times. In Section [3] 
we prove the compact embedding theorem. In the last section, we provide a new proof 
for Poincare-Friedrichs type inequalities for piecewise functions. Throughout the paper, C 
will be a generic constant that may depend on the domain Q and shape regularity /C of a 
triangle, of a triangulation, or of a family of triangulations. But otherwise, the constant is 
independent of triangulations. 



2. A TRACE THEOREM 

We first prove a trace theorem for piecewise functions. This result itself is a general- 
ization of a trace theorem of Sobolev space theory. It will be used in proving the compact 
embedding theorem, and in proving a Poincare-Friedrichs inequality by using the compact 
embedding theorem. We will need the following trace theorem on an element [2]. 
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Lemma 2.1. Let t be a triangle, and e one of its edges. Then there is a constant C depending 
on the shape regularity of t such that 



(2.1) 







-|e| j |Vn|2 


1 e 







Vug H\t). 



Theorem 2.2. Let Th he a shape regular, hut not necessarily quasi-uniform triangulation of 
Vt. There exists a constant C dependent on VL and the shape regularity constant of Th, hut 
otherwise independent of the triangulation such that 

(2.2) \\u\\L^i^dn)<C\\u\\Hi yueHl 

Proof. Let be a piecewise smooth vector field on Q whose normal component is continuous 
across any straight line segment, and such that cp-n = 1 on dQ. (The piecewise smoothness 
of (f) is not associated with the triangulation Th- A construction of such vector field is given 
below.) On each element t & Th, 'v/e have 

/ u^4>-n= div{u^(f)) = / (2mVm ■ + wMiv^). 

J dr J T J T 

Summing up over all elements of 7^, we get 



u 



an 



E 



QJe 



{2uVu ■ (f) + u^ div 0). 



If e is the border between the elements ti and T2 with outward normals rii and n2, then 
lu'^4>l = u\c()i-ni+u\cf)2'''^2i where Ui and U2 are restrictions of u on ri and r2, respectively. It 
is noted that although cf) may be discontinuous across e, it normal component is continuous, 



i.e., (f)^ ■ rii + (f}2 ■ 1^2 = 0. On the edge e, we have < 



0,oo,r2. 



Here, 



\ul - ul\ 



It is noted that = 2||u]{u}|, with {u} = {ui + U2)/2 being the average. We 



have 
(2.3) 



uw 



1/2 



nl/2 



< C'|0|o,oo,n 





1/2 






1/2 


jlnf 






-|er / |Vn|2 








Jse 


Jse 





Here, C only depends on the shape regularity of ti and T2. We used 5e = ri U T2 to denote 
the "co-boundary" of edge e, and we used the trace estimate (12. ip . It then follows from the 
Cauchy-Schwarz inequality that 



WWhidu) < C'(l<^lo,oo,n + I div(^|o,oo,n) 



\u\ 



+ 



\Vu\ 



Here the constant C only depends on the shape regularity of Th- The dependence on Q of 
the C in (12.21) is hidden in the (p in the above inequality. □ 
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We describe a construction of the vector field used in the proof. On the xy-plane, we 
consider a triangle OAB with the origin being its vertex O. Let the distance from O to the 
side AB be H . Then the field il^{x, y) = {x, y)/H has the property that i/^ ■ n = 1 on AB and 
V' ■ n = on OA and OB. Also \xl^\o,oo = max{|OA|, \OB\}/H and divV = 2/H. For each 
straight segment of dQ, we define a triangle with the straight segment being a side whose 
opposite vertex is in Q, then we define a vector field on this triangle as on the triangle OAB 
with AB being the straight side. We need to assure that all such triangles do not overlap. 
We then piece together all these vector fields and fill up the remaining part of the domain 
by a zero vector field. This defines the desired vector field used in the proof. 

3. Compact embedding in the space of piecewise functions 

We again assume that 7h is an arbitrary triangulation with shape regularity /C. For 6 > 0, 
we define a boundary strip of width 0(5). We let = f2 \ be the interior domain. 
The interior domain has the property that if a point is in Q^, then the disk centered at 
the point with radius 6 entirely lies in Q. We first show that when the strip is thin, the 
L'^{Qs) norm of the restriction on Qs of a function in Hj^ must be small. 

Lemma 3.1. There is a constant C depending on VL and the shape regularity /C ofTh, but 
otherwise independent ofTh, such that 

(3.1) / u\x)x<C6\\u\\li \/ueHI 

Jns 

Here Qs is a boundary strip of width 6 attached to dQ. 

Proof. We choose a piecewise smooth a vector field 4> whose normal component is continuous 
across any curve such that = on the inner boundary of Qs, and div (f) = 1 and \4>\ < C5 
on VLs. (A construction of such (p is given below.) We then extend by zero onto the entire 
domain VL. The extended, still denoted by 0, is a piecewise smooth vector field whose normal 
components is continuous over any curve in f2. We thus have 

/ = I div = ^ / wMiv = - ^ / 2m Vm ■ 4> + ^ ■ ^• 



The last term can be written as 



[ u^(f)-n= I M20-n+ V AmVI. 

Jdr Jan ^c-nJe 



Since the normal components of is continuous on edges in we use the same argument 
as in the proof of Theorem 12.21 cf . , (12.31) , to get 
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It then follows from Theorem 12.21 that 

(3.2) / < C|0|o,oo,Q 

The proof is complete since |0|o,oo,c < CS. The constant C depends on Q in terms of its 
interior angles and exterior angles at convex and concave vertexes, respectively. □ 

We describe a way to choose the boundary strip and construct the vector field (j) that 
was used in the proof. This field can be constructed by piecing together several special 
vector fields. We need some vector fields on rectangles, wedges, and circular disks. On the 
xy-p\a.ne, we consider the vertical rectangular strip R = (0,5) x (0,/). On this strip, we 
consider i/^^ = (x, 0). This vector field satisfies the condition that divi/>^ = 1, tl^j^ = on 
the left side, V/?, ■ 'n- = on the top and bottom sides and the maximum of is 6 that 
is attained on the right side. On a wedge W with its vertex at the origin, we consider the 
vector field ■j/'^y = {x,y)/2. This Vh/ satisfies the conditions that div il)^r = 1, t/^w ■ n = 
on the two sides of W, and IV'vkI ~ P/2 cit a point in W whose distant from the origin is 
p. On a circular disk C centered at the origin and of radius p, we consider the vector field 
il>Q = (1 — p^/r^)(x,?/)/2. Here r = (x^ + y'^Y^'^. This vector field satisfies the condition 
that div = 1 on the disk except at the center where it is singular. It points toward the 
center every where. And it is zero on the boundary. We use this field on a sector of the 
circle C, on the two radial sides of which we have i/j^; ■ n = 0. With these special vector 




Figure 1. Boundary strip near a convex vertex (left) and a concave vertex (right). 

fields, we can then assemble the on a boundary strip Qs- Along the interior side of each 
straight segment of dQ we choose a uniform strip of thickness 6. These strips overlap near 
vertexes. If Q is convex at a vertex, we introduce a wedge whose vertex is at the intersection 
of the interior boundary of the uniform strips, and whose sides are orthogonal to the meeting 
straight segments, see the left figure of Figure [1] If f2 is concave at a vertex, we resolve it by 
using a circular sector, centered near the vertex and outside of fl. The radius of the circle is 
slightly bigger than 6 such that the arc is continuously connected to the interior edges of the 
meeting strips, and the two radial sides are orthogonal to the meeting boundary segments. 
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see the right figure in Figure [T] With such treatment of the vertexes, the boundary strip 
is composed of rectangular strips attaching to major portions of straight segments of dQ, 
portion of wedges at convex vertexes, and portion of circular sectors at concave vertexes. 
See the shaded region in Figure [H We then transform i/^p,, ij^y^, and Vc various parts 
of VLs and assemble a </> on fi^. The vector field (p thus constructed is zero on the interior 
boundary of Vt^. Its normal components are continuous across any curve, and div(/) = 1 on 
VLs. The thickness of Vts is the constant 5 for the rectangular part. It is maximized to 5/ sin | 
at a convex vertex. It is minimized to a5 at the concave vertex, with < a < 1, a value can 
be chosen as, for example, 1/2. The norm \(f)\ has a maximum 5/ sin | at a convex vertex 
with 9 being the interior angle. Thus when 9 is small, |0| is big, and the estimate (13. 2p 
deteriorates. The norm |0| also has a local maximum near a concave vertex. It is bounded 
as 

, 1 - a sin I 

101 <c)- — ^. 

(1 — a) sm 2 

When the exterior angle is small this maximum is big. Also, one needs to choose a away 
from 1 and 0, to maintain a moderate thickness of the strip and a reasonable bound for |0| 
which affect the estimate (13. 2p . 

We then prove that functions in Hj^ are "shift-continuous" in L^, as stated in the following 
lemma. We extend a function u G H^^ to a function u on the whole by zero. 

Lemma 3.2. There is a constant C depending on VL and shape regularity fC of Th, hut 
otherwise independent ofTh such that 

(3.3) / [u{x + p) -u{x)]^x<C\p\\\u\\li MueHl 

Proof. Because for an element r G Th, smooth functions are dense in H^{t), we only need to 
prove (13. 3p for functions that are smooth on each element of Th- Let u be such a piecewise 
smooth function. Let p be an arbitrary short vector. We take 5 = |p| and choose a boundary 
strip Qs- The interior part of the domain has the property that if a; G then the line 
segment [x,x + p] G Q. We have 

[u{x + p) — u{x)]'^x = / [u{x + p) — u{x)]'^x + / [u{x + p) — u{x)]'^x. 



Using Lemma 13. we bound the second term as 



(3.4) / [u{x + p) - u{x)Yx < u\x)x<C\p\\\u\\j^i. 

We then focus on the first term. This integral can be taken on an equal measure subset of 
Q^. This subset is obtained by removing a zero measure subset that is composed of such 
point x: a; or a; + p is on an open edge e G S^, or the closed line segment [x, x + p] contains 
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any vertex of 7h, or [x, x + p] overlaps some edges of Sj^. By such exclusion, for any x in the 
remaining set, both the ends of [x, x + p] are in the interior of some open triangular elements, 
and [x, X + p] contains no vertex. The restriction of u on [x, x + p] is a piecewise smooth one 
dimensional function, which may have a finite number of jumping points in [x,x + p). By 
the fundamental theorem of calculus, we have 



u{x + p) — u{x) 



Vn(x + tp) ■ pt 



pe[x,x+p]n£'jl 

Note that the integrand in the integral may make no sense at t, if x + tp G S^, where u may 
jump. These points are excluded from the integration, where the jumping effect is resolved 
in the second term. On the segment [x, x + p], u may have a jump at p E [x,x + p] fl 
which is denoted by |n]p that is the value of u from the x side minus that from the x + p 
side. We thus have 

1 ^ 

_pe[x',x+p]n£'0 

When we take integral on (minus the aforementioned zero- measure subset), the first term 
is bounded as follows. 



[u{x + p) - u{x)Y < \p\ / |Vm(x + tp)|^t 

Jo 



no 



\P\ 



\Vu{x + tp)\hx = IpI 



|Vu(x + tp)pxt < IpI 



iVuPx. 



To estimate the jumping related second term, we write |u] 



|l/2|p 



-1/21 



U 



if p G 



[x, x + p]r\S^ is on the edge e, and use the Cauchy-Schwarz inequality to reach the following 
estimate. 

2 



u 



< 



en[x,x+p]=id 



^1 ^Men[i^,i^+p] 



en[x,x+p]j^d 



We show below that there is a C, depending on the domain Q and the shape regularity /C 
of Th , such that 



(3.5) 

We then have 



en[x,x+p]yt9 



no 



E 



u 



p€[x,x+p]n£'jl 



< c 



no 



E 

en[x,x+p]j^(!) 



|e| ^lufen[x,x+p]^- 



Every term in the right hand side is associated with a particular edge e E S^. Each edge 
e E is relevant to at most the points in the parallelogram fig in the Figure El Thus, by 
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changing the order of sum and integral, we have 




Figure 2. A p shift and fig for an edge e. 




Note that the second term may carry the smaller coefficient \p\ max{/i, |p|} such that the two 
terms are closer in the order. But we do not need such refined estimates. We thus proved 

(3.6) / [u{x + p)-u{x)]\<C\p\\\u\\jji \/ueHl 

The shift continuity (13.31) then follows from (13.61) and (13. 4p . We have shown that the set of 
zero extended functions is shift continuous in L^(M^). □ 

We give a proof for the estimate (13.51) . Let / be a straight line cutting through Q. Let Th be 
a shape regular triangulation with regularity constant /C. Then the sum of lengths of mesh 
line segments intersecting / is bounded independent of the triangulation. More specifically, 
we prove that there is a constant C, depending on the shape regularity /C, but otherwise 
independent of the triangulation Th such that 

(3.7) 1^1 ^ 

eefh and enl^0 
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We shall use some facts that follow from the shape regularity assumption. There is a min- 
imum angle 9]c among all angles of triangles of Th- The number of edges sharing a vertex 
is bounded by a constant C that only depends on /C. Let Ci and 62 be two edges shar- 
ing a vertex. There are constants Ci and C2 depending on /C such that |ei| < Ci|e2| and 
1 62 1 < C2|ei|. Without loss of generality, we assume / is horizontal. To simplify the argu- 
ment, we also assume that / does not pass any vertex. (This restriction can be removed by 
a slight modification of the following argument.) We first trim the set of intersecting edges 





I 




b / 


E 



























Figure 3. A line cutting through the triangulation. 



{e G £h and e fl / 7^ 0}. Consider the left most edge intersecting /, of which A is the end 
vertex shared by some other edges intersecting /. If A is above /, we examine all the edges 
intersecting / and sharing A in the counterclockwise order. We discard all such edges but 
the last one that is AB in Figure |3l (The next edge sharing A, as AC, does not intersect 
/.) The edge BC must intersect /. There could be other edges intersecting / and sharing 
the vertex B. We examine all the edges sharing B and intersecting / in the clockwise order. 
We discard all but the last one. (It is BC in the figure.) Now the vertex C is in the same 
situation as A, and we can determine the edge CD using the same rule as for AB. Then we 
determine DE, EF, and so forth. We do the trimming all the way to the right end of /. This 
procedure touches all the edges intersecting /, by either trimming an edge off or keeping it. 
The remaining edges constitute a continuous piecewise straight path as represented by the 
thick line in the figure. We denote this set by S^. It follows from the aforementioned facts 
about the shape regular triangulation that there is a constant C, depending on /C only, such 
that 

e€£^ and eDlj^d egf^ 

We consider a typical triangle bounded by / and S^, as the shaded one in the figure, whose 
sides are a, b, and c. Let the angle ZDEF be denoted by 6. We have 6 > die- Note 
that = + b"^ — 2abcos9. If 9 is obtuse, then a + b < \/2c. Otherwise, we have 



^2 _ t„2 



a- + 62)(1 _ cos^) + (a - bf cos9 < (a^ + b'^){l - cos9). Thus a + b < x/j^^c. In 
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any case, we have a + b < 



1— cos 8ic 



2 



c. We thus proved 



2 



1 — cos 6ic 



inn\. 



From this, (13. 7p follows. We can now prove the following compact embedding theorem. 

Theorem 3.3. Let be a (infinite) family of shape regular but not necessarily quasi- 
uniform triangulations of the polygonal domain Q, with a shape regularity constant /C. For 
each i, let Hi, be the space of piecewise functions, subordinated to the triangulation Th^, 
equipped with the norm (ll.ip . Let {ui} be a bounded sequence such that Ui G for each i. 
I.e., there is a constant C , such that \\Ui\\ui < C for all i. Then, the sequence {ui\ has a 



convergent subsequence in L'^{Q). 

Proof. It follows from (13. 3p that the sequence {ui} is a shift-continuous subset of L^(f2). The 
result then follows from the well-known condition for a subset of L^(f2) to be compact, see 



4. Poincare-Friedrichs type inequalities for piecewise functions 

Poincare-Friedrichs type inequalities have been generalized to spaces of piecewise func- 
tions. Several variants and proof methods of such inequalities can be found, for example, 
in [21 m [8]. We provide an alternative proof for this kind of inequalities by using the com- 
pactness result of the previous section, to demonstrate how to use the compact embedding 
theorem to obtain estimates that are uniformly valid with respect to triangulations. This 
method is a modification of the classical methods of proving some of the Poincare-Friedrichs 
inequality based on compactness argument. 

Theorem 4.1. Let Th be a shape regular, but not necessarily quasi-uniform triangulation of 
the polygon VL. Let f be a semi-norm defined on the space H\ that satisfies two conditions. 

1) There is a constant C that only depends on Vt and the shape regularity ofTh, otherwise 
it is independent of Th, such that 



2) For any constant c, f{c) = if and only if c = 0. 

Then, there exists a constant C depending only on the domain Q and the shape regularity 
constant ofTh, but otherwise independent ofTh, such that 



Theorem 2.12 in for example. 



□ 



(4.1) 




(4.2) 




< C 
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A major point here is that the constant C only depends on the shape regularity of Th- 
Had one only considered a particular triangulation, such inequality would easily follow from 
the Rellich-Kondrachov compact embedding theorem [10] and Peetre's lemma (Theorem 
2.1, page 18 in [9]). But this argument can not establish the independence of the C on the 
triangulation. To establish the independence of C of the triangulation, we need to consider 
the entire family of all triangulations that has a finite shape regularity constant as a family. 
We need to modify the arguments of proving the Peetre's lemma to handle the family of 
spaces. 



Proof of Theorem 4-1 ■ The inequality (14. 2 p is equivalent to 



(4.3) 



\u\ 



Hi 



< c 



[ iVnp + ^ / [uf + f\u) 



Vug Hi 



If there is not a constant only depending on the shape regularity of Th, such that this 
inequality holds, then there is a sequence of shape regular triangulations Th^ with a common 
shape regularity constant, and a sequence of functions Ui G Hi such that 



(4.4) 

and 

(4.5) 



\Ui\\Hl 



According to Theorem 13. 3[ this sequence of functions has a convergent subsequence, still 
denoted by Wj, in L'^iVL). We let the limit be Uq G L'^iVL). We claim that this Uq is a 
constant, and the constant is zero. We show this by verifying that the weak derivatives of 
Mo is zero. For a compactly supported smooth function 0, we have 



lim 



Uidicj). 



For an by using integration by parts on each triangle, we write the above right hand side 



as 



/ Uidi(j)= - ^ / <9lMi0+ ^ / [Mjni0- 



Here, on an edge e, = + u_ni_ if e is shared by r+ and r_ and ni+ is the first 

component of the unit outward normal of r+. Using Lemma [2. II to </>, it follows that 

1/2 



/ Uidicj) 


< c 


j 






Jn 
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In view of (14. 5p . we see that the weak derivative diUo is zero. For the same reason, ^2^0 = 0. 
Thus uq is a constant. We see from f l4.5p that — no||//i — ?■ 0. Thus, by the condition (14. ip 
f{uQ — Ui) < C\\Uq — Ui\\Hi — )■ as z — 7- 00. It follows from /(uq) < f{ui) + /(uq — Ui) V z 
that /(mq) = 0. Therefore, Mq = 0. Thus, UmiUhi ~^ when i ^ 00. This is contradict to 
f lOj) . ' □ 

In the classical Poincare-Friedrichs inequalities, typical forms of the semi- norm f{u) are 

1/2 



w 



f2{u) 



U 



or 



h{u) 



u 



Here F is, or a portion of, the boundary dVt, (which could be a segment with positive length 
of any curve in fi,) and u is, or a sub-domain of, Vt. Every one of these can be put in the 
position of the / in Theorem 14. II to obtain Poincare-Friedrichs inequalities for piecewise 
functions [21 H] . To see the validity of the inequality (14. 2 p for each of these semi- norms, one 
needs to verify the uniform boundedness condition (14. ip . in which the constant C is only 
allowed to depend on the shape regularity of Th and fi, and show that for any constant c, 
/(c) = if and only if c = 0. This latter condition is obviously met by all of them. For /i, 
the uniform boundedness condition follows from the trace theorem, see (12. 2p . For /2, the 
condition follows from the Holder inequality on F and the condition for /i. For /a, it follows 
from the Holder inequality that f?,{u) < |u;|^/^||m||l2(^) < 

By using a similar modification of the compactness argument, one can generalize the Korn 
inequality for shells (see |3]) to the space of piecewise functions, which seems useful for the 
analysis of discontinuous finite element methods for shells. 
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